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Abstract 

We study theoretically the behavior of vortices in a thin film superconductor placed close to a 
soft magnetic film. It is shown that the field from the vortex induces a magnetization distribution 
in the soft magnetic film, thus modifying the fields and vortex interactions. We suggest that the 
interaction between two otherwise identical vortices is attractive at short distances, but repulsive 
at larger distances. This is in contrast to the case without the soft magnetic film, where the force 
is always repulsive. 
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The structure of vortices in thin superconducting films was first investigated in detail by 
Pearl [1]. He found that such vortices interact mainly via their stray field, which extends 
far into the nonsuperconducting medium. Later, this theory has been extended to thin film 
systems with and without anisotropy ^ 0- | . More recently, the interaction between vortices 
and magnetic nanostructures has gained significant interest, due to the possible enhancement 
of the critical currents as well potential applications in future fluxtronic devices (see e.g. Refs. 

and references therein). So far the theoretical approaches developed 
have assumed hard nanomagnets, i.e. magnets that are not influenced by the vortex field 
(see e.g. Refs. and references therein). However, in many cases the magnetic film 

is softf7j, and the vortices will therefore induce a magnetization distribution in these films. 
The purpose of this paper is to make a first simple approach to study the vortex interactions 
in presence of a soft magnetic film that is homogenous in absence of vortices. 

To this end, we consider a simple model based on a thin superconducting film of infinite 
extent, located at z=0 with thickness d much smaller than the penetration depth of the 
superconductor. The surface is covered by a soft magnetic film of thickness smaller than that 
of the superconducting film. We also assume that there is no spin diffusion and proximity 
effects. In general, the current density is a sum of the supercurrents and magnetically 
induced currents, which can be expressed through the generalized London equation as^] 

V x J = ~H + -^V{p)6{z)e z + V x V x M , (1) 
A A 

where J is the current density, A is the penetration depth, H is the magnetic field and 
V x My is the magnetically induced current. Note that the magnetically induced currents 
are included as the last term on the right-hand side of Eq. ^ and are therefore generated 
in the same plane as the vortex. This is justified since we assume that the thicknesses of 
both the superconducting and magnetic films are smaller than the penetration depth, and 
the magnetic film is located very close to the superconductor. The vortex is aligned in the z 
direction, and its source function V(p) is assumed to be rotational symmetric. In the case 
of a Pearl vortex we may set V(p) = ($ // i o) ( Kp)> where $ is the flux quantum and fi the 
permeability of vacuum. 

Let us assume an initially homogenous soft magnetic film without domain walls, consist- 
ing of a single domain with in-plane magnetization in absence of external magnetic fields. 
Moreover, we assume that the free energy of this domain can be expressed as a sum of the 
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uniaxial anisotropy and the demagnetizing energy. The magnetic field from a vortex tilts 
the magnetization vector out of the plane due to the additional Zeeman energy, and it can 
be shown that for small tilt angles one has ||| 

M = (M p , M z ) « (m s , M s ^f] , (2) 



where M p = y/M* + M* and M z are the in-plane (radial) and perpendicular components 
of the magnetization, H p = «/ H% + Hy and H vz are the in-plane (radial) and perpendicular 
components of the vortex field, and M s is the saturation magnetization of the magnetic film. 
H a = M s — 2K U / hqM s is the socalled anisotropy field, where K u is the anisotropy constant 
of the magnetic film. Note that since we neglect the cubic anisotropy of the system, the 
in-plane magnetization direction must be the same as that induced by the vortex field. It 
is important to point out that the total magnetization vector M is not directed along the 
vortex field, since there is always a large component of the magnetization in the plane of the 
magnetic film due to the uniaxial anisotropy. Upon using Eq. we neglect the contribution 
from the exchange energy, which is justifiable for sufficiently large K u or small spatial field 
gradients. 

In order to solve the generalized London equation, we follow the method of Ref. fl. The 
current can only flow in a thin layer of thickness d (d <C A) , and it is therefore necessary to 
average over this thickness and consider the z component only 

/ Qj s 3J S \ 

H * + X e ( -nf - ) = V(p) + d\ e ( V x V x M) z , (3) 



dx dy 

where A e = X 2 /d is the effective penetration depth. The sheet current flowing in the thin 
layer is now given by J s = dJ. The Maxwell equation V X H = J gives 

j _dH ± _ dH^ j _ dH^ _ dHz 

dy dz ' y dz dx 

Since all derivatives d/dz are large compared to the tangential d/dp, we may set 

J,J. ~ H y — H y , jy ~ " , ^) 

where H^ and H^ (i=x, y) are the components at the upper and lower surfaces, respec- 
tively. Since the environments of the upper and lower half-spaces are identical, we have the 
boundary condition Hf = —H^, which results in 

J s x « -2H+, J s y « 2H+ . (6) 



Using V • H = 0, Eq. (jSJ) becomes 

377 

H z -2\ e -^ = V{p) + d\ e {V xV xM) z . (7) 

Using Eq. |2 an d the fact that our system has no volume charges, we find 
TT dH z fd 2 H z d 2 H z \ . , d\ e M s 

In order to solve Eq. (JHJ), it is useful to note that our system is rotational symmetric, and 
also that Vxff = V > H = outside the system. Therefore, we can introduce a scalar 
potential <fi which vanishes at z —>■ ±00 

I p 00 /* 00 

<Kp, z) = J J <j>{k) exp(ik ■ p - k\z\)d 2 k , (9) 

where k = ^fc 2 + fc 2 , and H = V</>. It is helpful to note that H z {k) = —k(f)(k) for the 
upper half-space. 

Applying the Fourier transform to Eq. (jHJ), we obtain 

^ = - (ITi^CTji ' (10) 

We assume that the vortex is located at the origin, and therefore the resulting scalar potential 



is 

"OO poo 



W J (27r) 2 7„ oo y_ 00 fc(l + 2A e A; - afc 2 ) 1 ' 

Due to the rotational symmetry the potential is found to be 



where we have used that 

-"•27T 



1>Z1T 

/ exp(ikpcos(J))d(J) = 2TrJ (kp) . (13) 
Jo 

To obtain the magnetic field components (in the radial and z direction), we apply the 
following formula 

^Mkp) = -kJx{kp) , (14) 



and find 
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= i[ kv ^ \ + uT- ak ^ - kWk ■ (16) 

In the case a = and V(k) = Q /u , the field is reduced to that of the standard Pearl 

fin 

solution |9j. We also notice that there is a divergency in k-space when 1 + 2X e k — ak 2 = 0. 
This is due to the fact that we assumed a very thin superconductor (d <C A). The divergency 
will therefore smoothen out upon solving the London equation for arbitrary thicknesses d, 
but this task is outside the scope of the current work. Nonetheless, it is seen that at small 
distances (large k) the magnetic field may change sign as compared to the standard Pearl 
solution. On the other hand, at large distances (small k) the field is basically not influenced 
by the soft magnetic film. A typical scale for the crossover is k ~ H a /M s d ~ 1/d, since in 
many practical cases M s ~ H a . We argue that the crossover follows this dependence also 
if we allow for a thicker magnetic film, but at some point one must take into account the 
magnetic volume charges. 

Let us now try to estimate the interaction energy and force between two vortices. To 
this end, we note that there are two contributions to this interaction. First, the energy 
associated with the two vortices (with indices 1 and 2, respectively) is given by 



E v =p d / H v V 2 d 2 p , (17) 

J —oo J —oo 

which for Pearl vortices is found to be 

EM = ^H z {p)d = ^ r k — J f p) dk . (18) 

27r/i Jo 1 + 2A e fc — ak 2 

One should also take into account the interaction between the vortex field and the vortex- 
induced magnetization, which can be found by using 



OO POO 



E m = -d / M r B 2 d 2 p . (19) 



oo J — oo 



Note that the magnetic induction in the magnetic film is B 2 = po{H 2 + M 2 ), which means 
that to the first order E m is a constant proportional to poM 2 . It should be emphasized that 
this is only correct when the vortex field is substantially weaker than the anisotropy field 
and magnetization of the magnetic film. 

Based on the above observations we estimate the force between two vortices to be 

yF ' 27cpoJ l + 2\ e k-ak 2 v ; 
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When the distance is large (p ^> A e ), the soft magnetic film does not influence the vortex in- 
teraction, and the force is taefore^ive and govern*, by ft e standard Pear, Solutlon Q. 
On the other hand, at smaller distances (still larger than the coherence length £), we may 
approximate the force by 

F( P )~-^,(< P «^d, (21) 
2np ap H a 

where we have used that J °° J\{kp)dk = 1/p. Therefore, at small distances p the force is 
attractive in presence of a soft magnetic film, which is quite surprising. However, one can 
interpret this as a result of the currents generated by the vortex-induced magnetization. 
These currents generate a magnetic field which opposes that produced by the vortex in 
absence of a magnetic film, and this field interacts with the core of the second vortex. For 
comparison, we note that in absence of a soft magnetic film the force at small distances is 
given by 

F(p) « , 2 , i < P « A e , (22) 

which is seen to be repulsive with a magnitude that is decaying faster with increasing p than 
Eq. EU 

One may expect that the presence of a soft magnetic film results in new vortex con- 
figurations and also influences the creation of vortices at the superconducting transition 
temperature. We therefore hope that this study will stimulate further work in this field. 
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